Controlling spin current in a trapped Fermi gas 
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We study fundamental features of spin current in a very weakly interacting Fermi gas of 6 Li. By 
creating a spin current and then reversing its flow, we demonstrate control of the spin current. This 
reversal is predicted by a spin vector evolution equation in energy representation, which shows how 
the spin and energy of individual atoms become correlated in the nearly undamped regime of the 
experiments. The theory provides a simple physical description of the spin current and explains 
both the large amplitude and the slow temporal evolution of the data. Our results have applications 
in studying and controlling fundamental spin interactions and spin currents in ultra-cold gases. 



Spin dynamics and spin currents have been extensively 
studied in condensed matter physics [l[ . Active manipu- 
lation of the electron spin can be used for data process- 
ing and storage For example, a spin current can be 
used to excite or reverse the magnetization of a nanomag- 
net 0,0, HI- In ultracold atomic physics, spin-current- 
related phenomena have been observed both in a Bose 
gas Q and in a Fermi gas Q. We report on the origin 
and control of spin current in a weakly interacting Fermi 
gas. 

An optically trapped Fermi gas of 6 Li is a rich system 
in which the strength of interactions between atoms can 
be controlled by applying a variable bias magnetic field 
tuned near a Feshbach resonance [a, [9(. Close to res- 
onance, the Fermi gas exhibits strong interactions fioj ] . 



which have been widely studied [ll| . This regime of- 
fers unprecedented opportunities to test nonperturbative 
quantum many-body theories. In contrast, little investi- 
gation has been done for a weakly interacting Fermi gas, 
where the s-wave scattering length, ai2 , for atoms in op- 
posite spin states can be tuned smoothly from small and 
positive to small and negative. In this regime, spin segre- 
gation is observed, where atoms of one spin move outward 
in the trap, while atoms with the opposite spin move in- 
ward. Previous observations of spin segregation Q have 
shed new light on the study of this regime. 

An overdamped spin wave theory 12, 13, has been 
used to explain the spin segregation observed in a Bose 
gas of 87 Rb confined in a magnetic trap [f|. In those 
experiments, the collision rate between the atoms was 
large compared to the axial trap frequency and sufficient 
to ensure a thermal momentum distribution. The pre- 
dictions are in good agreement with the measurements 
for the Bose gas. In contrast, the corresponding theory 
of overdamped spin waves for a Fermi gas disagrees with 
the experiments by two orders of magnitude in ampli- 
tude, and predicts an oscillation of the density profile, 
which is not observed Q- 

In the Fermi gas experiments with 6 Li, collisions be- 
tween atoms in the same state are prohibited due to the 
Pauli principle, and the scattering length between atoms 
in opposite spin states is magnetically tuned to be very 
small. In this case, velocity changing collisions between 



atoms in different states occur at a rate of only ~0.3 Hz, 
small compared to both the axial trap frequency and the 
spin segregation rate. We explain the observed spin seg- 
regation in this regime by a nearly undamped spin wave, 
in which the spin vector of each atom is correlated with 
its energy. 

To understand the origin of the spin-segregation and 
corresponding current, consider the evolution of the spin 
vectors of atoms vibrating almost freely along the axial 
(long) direction of a cigar-shaped optical trap. Atoms in 
states |t) and ||) oscillate in the trap with frequencies ui x ^ 
andw x i, respectively. For the trap, ZJ X = (uj x ^+uj x i)/2 = 
2n x 145 Hz. As the magnetic moments of the two spin 
states are not identical, the finite curvature of the bias 
magnetic field causes a small difference in the axial con- 
fining potentials and hence in the oscillation frequencies 
for the two spin states, Slo x = lj x i — uj x -\ = —2ir x 2.5 
mHz (the corresponding difference in the transverse os- 
cillation frequencies in the optical trap is negligible, due 
to the tight transverse confinement). The small differ- 
ence in the axial frequencies correlates the precession rate 
Q(E) of an atomic spin vector in the x — y plane with 
the energy of the atom. This along with binary collisions 
causes a spin wave. 

Fig- Ha) shows that the resonance frequency tv res for a 
radio-frequency (rf) transition between states n, | ) and 
|n, |) is shifted by (n + 1/2) (w^ — u> x i), where n is the 
harmonic oscillator quantum number for the axial direc- 
tion, which does not change in the transition. The axial 
energy E of the atom determines n + 1/2 = E/fuD x . As 
the collision rate for the Fermi gas of 6 Li in the experi- 
ment is very low compared to the axial trap frequency, 
collisions do not significantly change the energy of each 
atom over the time scale of the segregation, and the pre- 
cession rate can be written as 



n{E) = -{5lu x /uj x )E/?i. 



(1) 



As a result of this energy-dependent precession rate, the 
magnitude of the precession angle of the spin vector in 
the x — y plane is larger for atoms with high energy than 
for atoms with low energy. 

When two coherently-prepared atoms collide, the 
energy-dependent precession angle then leads to a corre- 
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FIG. 1: Spin-wave formation, (a) For an rf transition be- 
tween harmonic oscillator-spin states | n, j) and | n, J,), the 
resonance frequency u res (denoted by the red arrow) increases 
with n, due to the difference in the harmonic oscillator fre- 
quencies for the two spin states. An rf pulse initially cre- 
ates x-polarized spins (in the rotating frame), (b) The spin 
vector for atoms of high energy E > precesses more than for 
atoms of low energy E < . Binary collisions then cause the spin 
vectors to rotate about the total spin vector S, producing a 
z-polarized spin wave. 



lation between the z— component of the spin vector and 
the energy. The collisional interaction results in a ro- 
tation of each atom's spin vector about the total spin 
vector, which is conserved. The sense of the rotation is 
determined by the sign of the scattering length ai2, and 
the relative angle of the spin vectors. As both atoms have 
spins in the x-y plane, the rotation of each spin about the 
total spin in the x-y plane produces spin components out 
of the x-y plane, as shown in Fig. W[b). For a positive 
(negative) scattering length, atoms with higher energy 
E > will accumulate a negative (positive) z— component, 
while atoms with lower energy E < will accumulate a posi- 
tive (negative) z— component. This process correlates the 
z— component of the spin with the energy, i.e., S Z (E). 

The spin density vector in coordinate space is then 
determined by the axial harmonic oscillator wavefunc- 
tions 0_e(x) as S(x,t) = J dE S(E, t)\(j>E(x)\ 2 , where we 



assume that there is no coherence between different en- 
ergy states. The energy-dependent spin vector S(E,t) 

is determined from the Heisenberg equations S(E,t) — 
(i/h)[H, S], where the components of S(E, t) are written 
in terms of creation and annihilation operators in energy 
representation. 

In a one-dimensional approximation, the Hamiltonian 
operator for a Fermi gas in the optical trap is 



H — Hq + Hi, 



(2) 



In a frame rotating at the unshiftcd hypcrfmc transition 
frequency, 

H = Y,e(n 1 (E)+N i (E)) +J2mE)S z (E), (3) 



where N^i(E) are the number operators for each state 
a ndS,(E) = [N^E)-N l (E)}/2. 

Collisions produce a contact interaction between op- 
posite spin states. Averaging over the transverse coordi- 
nates {z, y), the collision operator takes the form 



Hi, 



dx^Ux^Ux^iix^ix). (4) 



Here, m is the atomic mass and a p is radial 1/e width 
for a fit of a Gaussian distribution of the trapped cloud. 
Eq. 2] can be written in energy representation using 
i(>1 t l(x) = J dE a-f ; i(E) 4>e(x), where 4>e(x) is the ax- 
ial harmonic oscillator wavefunction. 

Using S(E,t) = (S(E,t)), we obtain the evolution 
equations 



dS(E,t) 
dt 



n{E) x S(E,t) 

J dE 1 g(E', E) S(E', t) x S(E, t), (5) 



where 



l(E',E) = - 
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(6) 



Here E m i n = min(E, E') and we have assumed a WKB 
approximation [l5| for (J>e{x) between the classical turn- 
ing points, since the energies are in the classical regime. 

Eq. [5] is a primary result of this paper. The first term 
describes the energy-dependent spin precession, while the 
second term describes the rotation of the spin vector 
arising from binary collisions. The initial spin vector is 
S x (E,t = 0) = Nexp(-E/K B T)/(2k B T), where N is 
the total number of atoms and the factor 1/2 arises from 
the definition of the spin vectors. Fig. [5] shows that the 
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FIG. 2: (a) Difference between the spin-down and spin-up 
densities nj,— n-j at 220 ms in units of no = (n^o+ n To)/2- Here 
riio is the initial spin density for each state at the trap center 
before spin segregation occurs. Data (blue line); Theory (red 
curve). The data is taken for di2 = — 4.5ao; (b) nj —n-f at the 
trap center versus time for ai2 = — 4.5ao. Data (blue dots); 
Theory (red curve). 



predictions obtained by numerical integration of Eq.[5]are 
in good agreement with the data and trap parameters of 
Ref. 0. 

Our description of the spin vector evolution in energy 
representation can be compared to two recent theories 
based on a collisionless Boltzmann equation [l6|, 17 1. 
These approaches provide a phase space description in 
one dimension for a weakly interacting two-component 
Fermi gas. Their results are in very good agreement with 
the predictions of Eq. O both in amplitude and temporal 
evolution. 

The spin-energy correlation description is particularly 
useful when the evolution is nearly Hamiltonian, where 
it provides a very simple physical picture for a weakly 
interacting spin system. Eq. makes it apparent that 
the spin current can be reversed, causing the spatial dis- 
tribution to return to the unsegregated state. Two steps 
are required: A it pulse is applied to reverse the sign of 
S z and either S y or S x ; the sign of the scattering length is 
inverted by sweeping the bias magnetic field through the 
zero crossing. Then the cloud will start to merge. The 
quantitative predictions of our numerical simulations are 



shown in Fig. [3ja). Our experiments confirm this pre- 
diction. Immediately after the rf it pulse and inversion 
of the sign of the scattering length, the difference in the 
spin-up and spin-down densities at the cloud center is 
30% of the average density. As the cloud merges, this 
difference decreases to zero in ~ 70 ms, in contrast to 
the case without current reversal, where spin segregation 
persists for a few seconds. Then the density difference 
continues to evolve and returns to 30%. 

To perform these experiments, we prepare a spin seg- 
regated sample of 6 Li Fermi gas at a ±2 — 8.1ao [18| (the 
bias magnetic field B = 529.8 G is calibrated with rf 
spectroscopic techniques). A sample of 6 Li atoms in a 
50-50 mixture of the two lowest hyperfine states is loaded 
into a CO2 laser trap with a bias magnetic field of 840 
G, where the two states are strongly interacting. Evap- 
orative cooling is performed to lower the temperature of 
the sample [10(. The magnetic field is then increased in 
0.8 seconds to a weakly interacting regime at 1200 Gauss 
where an on-resonance optical pulse of 40 /is is applied 
to remove atoms of one state, while leaving atoms in the 
other state. With a single state present, the magnetic 
field is lowered in 0.8 seconds to 529.8 G. Then a 40 ms 
rf pulse (center frequency 75.613 MHz and sweep range 
35 kHz) is applied on the |t}— ||) transition to create a 
50-50 coherent superposition of the two spin states. Note 
that, as the frequency passes through resonance, coher- 
ence is created on a time scale of a few milliseconds, short 
compared to the time for the total sweep and for spin 
segregation to occur. At the final optical trap depth, the 
measured trap oscillation frequency in the transverse di- 
rections is oj± = 2tt x 3900 Hz, while the axial frequency 
is uj x = 2ir x 120 Hz. The total number of atoms is 
N ~ 4.0 x 10 5 . The corresponding Fermi temperature is 
Tf — 6 /iK. The sample temperature is T ~ 33 /iK. The 
peak atomic density at T is 6 x 10 11 /cm 3 . The axial and 
radial 1/e widths for a fit of a gaussian distribution to 
the initial density profile of the sample are ~ 400 /im and 
~ 12 /xm, respectively. 

At 40 ms after the first rf pulse, when difference in 
the spin-up and spin-down densities is 30% of the av- 
erage density, we change the bias magnetic field from 
529.8 G (012 = 8.1a ) to 525.2 G (a 12 = -8.1a ) in 5 
ms. Then we apply an rf it pulse (duration 40 ms, center 
frequency 75.596 MHz and sweep range 4 kHz), which 
flips the spins, as shown in the 0~ and + images of 
Fig. [3jb). Finally, we take absorption images of atoms 
in both states (in separate experimental cycles) at var- 
ious times after the rf it pulse, Fig. [3] [19(. The entire 
time sequence is done in random order and repeated 6 
times. The error bars are statistical and arise from run- 
to-run variations in the atom number, magnetic field and 
excitation frequency. 

As a further test of this idea, we applied either the 
scattering length sign change or the rf it pulse, but not 
both, and looked for merging followed by segregation. 
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FIG. 3: Reversal of the spin current in an ultracold Fermi 
gas. (a) Difference in the spin densities nj, — n-f at the trap 
center versus time, showing return to at t ~ 70 ms, after 
an rf 7r pulse is applied and the sign of scattering length is 
inverted at t ~ 0. Prediction (red line); Data (blue dots), 
(b) Images of the spin waves corresponding to the time in 
(a). Each image is the z-component of the spin density as a 
function of position x along the axial direction of the trap. 
0~ ms corresponds to time just before the rf 7r pulse and + 
ms corresponds to time just after the rf ir pulse and reversal 
of the scattering length. 

We found no reversal of spin segregation. This verifies 
that both operations are required to observe the reversal 
of spin segregation, as predicted. 

In an additional experiment, we increased the length of 
the initial spin segregation time, to determine the longest 
time scale r over which the spin segregation is reversible. 
We expect that r must be significantly smaller than the 
velocity changing collision time (~ 3 sec) for reversal to 
occur and find r = 200ms. 

Our experiments suggest that broad manipulation 



of the spin dynamics and creation and study of non- 
equilibrium systems in arbitrary spin mixtures is possi- 
ble, by using general rf pulse sequences and by temporally 
varying the scattering length and the spatial dependence 
of the magnetic field. 
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